1 



Relaxation times distributions for an anomalously diffusing particle 

Noelle POTTIER* 

Laboratoire Matiere et Systemes Complexes (MSG) 

O : UMR 7057, CNRS & Universite Paris Diderot 

(N ■ 

F-75205 Paris Cedex 13, France 



in 

(N 



in 



X 



Abstract 



As well-known, the generalized Langevin equation with a memory kernel decreas- 
ing at large times as an inverse power law of time describes the motion of an 
! anomalously diffusing particle. Here, we focus the interest on some new aspects 

! of the dynamics, successively considering the memory kernel, the particle's mean 

• velocity, and the scattering function. All these quantities are studied from a unique 

^ ■ angle, namely, the discussion of the possible existence of a distribution of relaxation 

^ . times characterizing their time decay. Although a very popular concept, a relax- 

^ ! ation times distribution cannot be associated with any time-decreasing quantity 

O I (from a mathematical point of view, the decay has to be described by a completely 

monotonic function). 

Technically, we use a memory kernel decaying as a Mittag-LefHer function (the 
>• . Mittag-LefHer functions interpolate between a stretched or compressed exponen- 

tial behaviour at short times and an inverse power law behaviour at large times). 
I We show that, in the case of a subdiffusive motion, relaxation times distributions 

■ can be defined for the memory kernel and for the scattering function, but not for 

^ ■ the particle's mean velocity. The situation is opposite in the superdiffusive case. 
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1. Introduction 



In this paper, we study some new aspects of the dynamics of an anomalously 
diffusing particle whose motion is described by the generalized Langevin equation: 



where v{t) is the particle velocity, F{t) is the Langevin force acting on the particle, 
as modelized by a stationary Gaussian random process of zero mean, and 7(t) is 
a dissipative memory kernel. Only the knowledge for t > of this latter function 
is required in Eq. (1.1), but, when necessary, 7(t) will be defined for any t as an 
even function of time. As well-known, for a particle evolving in a bath at thermal 
equilibrium at a temperature T, the consistency of the above model requires that 
F{t) and j{t) are not independent quantities, but instead are linked by the second 
fluctuation-dissipation theorem, that is: 



where ks is the Boltzmann constant. (The symbol (. . .) denotes the average over 
the realizations of the noise). 

In a previous paper [1], we studied this problem in detail in the case of a 
retarded memory kernel defined by its algebraic decay at large times {■^{t > 0) oc 
^ with 0<5<lorl<5<2). We showed in particular that the mean-squared 
displacement of the particle then varies at large times as: 



the motion of the particle being thus subdiffusive in the case < 6 < 1, and 
superdiffusive in the case 1 < 5 < 2. The subdiffusive case is especially interesting 
in relation with experiments in condensed matter physics. For instance, the above 
described model has been used to analyse the subdiffusion process undergone by 
the distance between a donor and an acceptor of electron transfer within a single 
protein molecule [2]. 

In the present paper, following [3], we will use a more general model for the 
retarded memory kernel, namely, a Mittag-Leffier function of index 5 with the 
same large time algebraic decay as in [1] and a modified exponential behaviour 





{F{t)F{t'))^mkBT^{t-t'), 



(1.2) 



(Ax2(^))~^^ 



(1.3) 
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at short times. In this framework, we will focus the interest on some new aspects 
of the dynamics of the anomalously diffusing particle, successively considering the 
memory kernel, the particle's mean velocity, and, last, the scattering function. All 
these quantities will be studied from a unique angle, namely, the discussion of the 
possible existence of a distribution of relaxation times characterizing their time 
decay. The concept of a relaxation times distribution is indeed a very popular 
one, but the quantity of interest (i.e., for which the existence of a distribution 
of relaxation times is discussed) needs to be made precise. For instance, it is not 
obvious that, in a given range of values of S, such a relaxation times distribution 
exists for each one of the above quoted physical quantities. 

The paper is organized as follows. In Section 2, we recall the general con- 
dition for the existence of a distribution of relaxation times associated with any 
time-decaying function, and the way this distribution can be obtained. In Sec- 
tion 3, we focus the interest on the retarded memory kernel 7(t) as described by 
a Mittag-Leffler function of index d. We propose a microscopic model justifying 
such a choice. To this end, we show how, in the spirit of the Mori-Zwanzig mem- 
ory function theory, the Mittag-Leffler memory kernel 7(t) can itself be viewed 
as evolving according to a Langevin-like equation associated with a power law 
second-order memory kernel, a feature which, to the best of our knowledge, does 
not seem to have been noticed before. The possible existence of a relaxation times 
distribution associated with ^{t) is then discussed in detail. The particular forms 
that this distribution (when it exists) takes in the short and long time regimes 
are discussed. In Section 4, we turn our attention to the particle's mean velocity 
{v{t)). We recall its general expression as found in [3], as well as its approxi- 
mate form in the long time regime as given in [1]. Then we use this latter ex- 
pression to discuss the possible existence of a relaxation times distribution as- 
sociated with {v{t)). In Section 5, we study the dynamics of the particle over 
a coarse-grained time scale corresponding to the description of the motion by 
a non-inertial Langevin equation. We introduce the scattering function Fs{q,t), 
defined as the spatial Fourier transform of the diffusion front, and we discuss 
the possible existence of a related relaxation times distribution. In Section 6, we 
provide an interpretation for our results, and we summarize our conclusions. 
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2. General condition for the existence of a distribution of relax- 
ation times 

In some instances, the time decay of the physical quantities of interest is described 
by completely monotonic functions. A function f{t) is said to be completely mono- 
tonic in the interval (a, b) if we have [4]: 

f^''\t)>0, a<t<b, A; = 0,1,2... (2.1) 

We will assume f{t = 0) = 1 in what follows. 

To begin with, let us recall the Bernstein's theorem, according to which the 
necessary and sufficient condition for a function f{t) to be completely monotonic 
in the interval < t < oo is that it can be written as [4] : 

POO 

f{t)^ / e-^'da{k), (2.2) 

where a{k) is bounded and non-decreasing, and the above integral converges for 
< t < oo. If a{k) is a function with a continuous first derivative /u(/c), f{t) appears 
as the direct Laplace transform of ii{k) > 0. The function //(/c), normalized such 
that Jg°° iJL{k)dk = 1, can be viewed as a probability density of rate constants [5]. 

2.1. Distribution of relaxation times 

From Eq. (2.2), we can obtain the distribution of relaxation times associated with 
the time decay of the completely monotonic function f{t). Setting: 

/>oo 

f{t)= / P(r)e-*/^dT, (2.3) 
we get from Eq. (2.2) and the definition fj,{k) = da{k)/dk: 

P(r) = ^^{l} (2.4) 

The function P(t), normalized such that P{T)dT = 1, is a distribution (more 
precisely, a probability density) of relaxation times. Conversely, we have: 
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2.2. Getting P{t) from the Laplace transform of f{t) 

In some cases, it may be more convenient to start, not from f{t), but from its 
Laplace transform f{z), as defined by: 

/>)=/ f{t)e-^'dt. (2.6) 

The decaying function f{t) can be obtained from f{z) by inverse Laplace trans- 
formation: 

/W = ^^/>K*c/z, (2.7) 

where Br denotes the Bromwich path in the complex 2;-plane, i.e. a line "Rez = a 
with (7 larger than the abscissas of the possibly existing singularities of J{z) and 
'ismz running from — oo to +oo. For those functions f{z) that have a cut along 
the negative real axis, it is possible to bend the Bromwich path into a Hankel 
path, i.e. a loop which starts from — oo along the lower side of the negative real 
axis, encircles the origin in the positive sense, and ends at — oo along the upper 
side of the negative real axis. If f{z) has poles in the first Riemann sheet, their 
contribution has to be taken into account, and it must be added to the contribution 
of the Hankel path. 

This latter contribution to f{t) reads: 

~ $>m|/(a;e*'^)}e-^*da;. (2.8) 

Thus, when no poles are present, the function: 

-^^m[f{ke'^)] (2.9) 

can be identified with the probability distribution of rate constants: 

^,{k) = ~'~Sm[f{ke^'')], (2.10) 

from which P{t) can be deduced using Eq. (2.4). 

3. The Mittag-Leffler retarded memory kernel 

Let us now come to the above described purpose, that is, the discussion of some 
relaxation times distributions intervening in anomalous diffusion. Following [3] , we 
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will introduce in the generalized Langevin equation (1.1) a Mittag-Leffler retarded 
memory kernel, namely: 

5-, 



lit) = ^Es 







-ii) 



t > 0, 



(3.1) 



with < 5 < 2. In Eq. (3.1), to acts as a characteristic memory time linked to 
the particle's environment (thermal bath), and 75 is a proportionality coefficient 
depending on the exponent 5 but independent of time. The Ea{y) function denotes 
the Mittag-Leffier function [6]-[8] defined through the series: 



Ea{y) = J2 ^, \ , a>0, 
^ T{an + 1) 

n=U 



(3.2) 



where T denotes the Euler Gamma function. For 5 = 1, we just have an ordinary 
exponential decay as pictured by j{t) = (7i/to) exp(— t/to)- 

We assume 5 7^ 1 in the rest of this subsection. As underlined for instance in 
[9], the choice (3.1) is especially convenient, since it allows ^{t) to interpolate, from 
a "stretched exponential" (in the case < 5 < 1) or a "compressed exponential" 
(in the case 1 < (5 < 2) for t <^ to, to an inverse power law for t ^ to. More 
precisely, retaining only the two first terms of the series expansion (3.2), we can 
approximately write, in the short time regime [9], [10]: 



-J exp 



(t/to) 



T{S + 1) 



t < to. 



(3.3) 



In the long time regime, we have [6] : 

^^^^ I (t^) 



-5 



t > to. 



(3.4) 



Note that, in this latter regime, we recover the power law retarded memory kernel 
previously used in [1]. Interestingly, the characteristic memory time to linked to 
the bath is in fact not involved in the long time expression for 7(t) (Eq. (3.4)). 
Setting 75 = ujg~\ Eq. (3.4) reads: 



lit) 



00} 



T{l-6) 



{u>sty 



t > to- 



(3.5) 



Formula (3.5) shows that the pertinent time scale is then Ug ^, which for consis- 
tency we assume to be much larger than to. Summarizing, two different time scales 
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enter into play in the Mittag-Leffler memory kernel (3.1), a short one denoted by 
to, and a long one denoted by ojg^- 

3.1. The short time regime 

Let us discuss in more detail the short time regime of the Mittag-Leffier function 
Es[—{t/to)^], as given by the modified exponential exp[— in formula (3.3). 

As underlined in [5] , the denominations of "stretched" and "compressed" ex- 
ponentials are somewhat of a misnomer in the short time regime which precisely 
interests us now. Indeed the so-called stretched exponential exp[— in the 
case < 5 < 1 is actually a function whose decay with time is faster than that of 
the ordinary exponential exp[— r(i+i) ] < t < to (but slower afterwards, which 
is in fact the reason for the denomination of "stretched" exponential). The initial 
decay rate of the retarded memory kernel (3.3), which is given by the derivative 
[d'y{t)/dt]{t = 0), is even equal to — oo in this case. 

In the case 1 < 5 < 2, the time decay of the so-called compressed exponen- 
tial exp[~ i^^s+i-j ] is slower than that of the ordinary exponential exp[— p^^^^] for 
< t < to (but faster afterwards, hence the denomination of "compressed" ex- 
ponential). The initial decay rate of the retarded memory kernel vanishes in this 
case. 

3.2. The Laplace transform of the Mittag-Leffler memory kernel 

Following the same procedure as in [1], we will solve the generalized Langcvin 
equation (1.1) using Laplace transformation. Introducing the Laplace transform 
^{z) of the retarded memory kernel ^{t) defined in the usual way, namely: 




(3.6) 



we get from Eq. (3.1) (as shown for instance in [8]): 



7(2;) 



1 



(3.7) 



z + tQ z 



For \z\<t:t. 



-1 
' 



we have: 



7(2;) ~ 752;' 



.5-1 



(3.8) 



that is: 



7(2) col'z'-' 



(3.9) 
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an expression consistent with the choice made in [1] , as well as with the large time 
behaviour of ^{t) as given by Eq. (3.5). 

3.3. A microscopic model for the Mittag-LefHer memory kernel 

In the spirit of the Mori-Zwanzig memory function theory [11]-[12], let us assume 
that the time evolution of the Langevin force F{t) acting on the particle velocity 
(a force which is denoted by Fi (t) for more clarity in what follows) can in turn be 
described by the Langevin-like equation: 

A{t)+ f 4>{t-t')F^{t')dt' = F2{t), (3.10) 

where (j){t) is a second-order memory kernel associated with the time evolution of 
F\ (t) , and F2 (t) acts as a second-order random "force" . Assume in addition that 
F2 {t > 0) is uncorrelated with Fi (0) . The correlation function Ci (t) = {Fi {t)Fi (0)) 
then obeys the homogeneous evolution equation: 

Ci{t)+ [ (l){t-t')Ci{t')dt' = 0, t>0. (3.11) 
Jo 

The memory kernel ^{t) = Ci (t) /mkBT of the generalized Langevin equation (1.1) 
obeys the same evolution equation: 

j{t)+ [ 4){t-t')^{t')dt' = 0, t>0. (3.12) 
Jo 

We want to solve Eq. (3.12) for a given ^{t = 0). Applying Laplace transformation, 
we get: 

^7(2) -7(t = 0) + 0(^)7(z) =0. (3.13) 
We thus recover for 7(z) an expression similar to formula (3.7): 

7(^) = (3.14) 

z + (p{z) 

More details on the derivation of Eq. (3.14) can be found in Appendix B. Compar- 
ing formulas (3.7) and (3.14), makes it clear that the choice of the Mittag-Leffler 
memory kernel ^{t) as given by Eq. (3.1) corresponds to the modelization: 



<^(z) = ^oV-^ 



(3.15) 
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that is, to an inverse power law second-order memory kernel (f){t): 

Since 5 < 2, the function decreases algebraically towards zero at large times 
{t » to), except for 5 = 1, in which case this algebraic tail vanishes, leaving us 
with the delta function behaviour (f){t) = {l/to)S{t), in accordance with the formal 
representation of the Dirac generalized function 5{t) = t~^/T{0), t > [13]. Note 
that only the short time scale to is actually involved in the expression for (but 
not the long time scale oj^^). We have (/)(t) < for < 5 < 1, and (/)(t) > for 
1<S<2. 

3.4. The relaxation times distribution associated with ^{t) 
The question of the existence of a relaxation times distribution associated with 
a Mittag-LefHer decay has been extensively discussed in the literature [10], [14]- 
[15]. We only summarize here the main results. Following [8], we write 7(t) as the 
inverse Laplace transform of j{z) as given by formula (3.7): 

7W = |^/^^r^e-d.. (3.17) 



where Br denotes the Bromwich path in the complex 2-plane {^ez = a with 
c ^ ^0^^)- The Bromwich path is then bent into the equivalent Hankel path Ha(£), 
i.e. a loop which starts from — oo along the lower side of the negative real axis, 
encircles the origin in the positive sense along a circle of small radius £, and ends 
at — oo along the upper side of the negative real axis. In both cases < 5 < 1 and 
1 < 5 < 2, we obtain [8]: 

7(^) = 7i(i) + 72(t), t>0, (3.18) 
where 7i(t) is the contribution of the Hankel path: 

7.W = |^/ ^r^e-d., (3.19) 



Ha(e) Z + tQ Z 

and 72 (t) is the contribution of the relevant poles, i.e. the poles situated in the 
first Riemann sheet, of l/(z + tQ^z^~^). 
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In the case < 5 < 1, no such poles exist [8], and ^{t) = 7i(t) relaxes towards 
zero without oscillating. The result takes a form analogous to Eq. (2.2) (except for 
the fact that we do not have ^{t = 0) = 1, but instead ^{t = 0) = ^s/to)'- 

/•OO 

7(t) = ^/ i^^{k)e-^'dk, (3.20) 
with the rate constants distribution fi-yik) as given by: 

(3.21) 

A distribution of relaxation times can thus be defined in this case (these relaxation 
times being linked to the thermal bath), which means that we can write: 

7(0 = ?/ P,{r)e-'/-dT, (3.22) 
f^o Jo 

with P-yir) a non- negative function normalized such that P^{T)dr = 1. From 
Eqs (2.4) and (3.21), we get [8],[10],[13]-[14]: 

P,(r) = ^ i , (3.23) 

\to/ \tQ/ 

The distribution Pj{t) as given by Eq. (3.23) is non-negative and normalized: 

POO 

/ P^(t) dr = Es{t = 0) = 1. (3.24) 
Jo 

It is widely referred to in the literature as the Cole- Cole distribution of relaxation 
times [16]. Its short time and large time limiting forms are discussed in detail in 
Appendix A. The behaviour of P-yij) for r <^tQ is especially interesting. We have 
an integrable divergence of the distribution of small relaxation times: 

_ , , sin^TT 1 / T \ 5-1 

P^(^) ^ __ _ j , r^to. (3.25) 

In the case 1 < 5 < 2, the situation is more intricate since ^{t) contains in 
addition an oscillating part. It can be written as [8]: 

^W = i +jy.(r)e-"^^r], (3.26) 
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where the function P'yir), which is still given by Eq. (3.23), cannot be given the 
meaning of a relaxation times distribution. We now have: 

r°° 2 2 

P^{T)dT = Es{t = 0)--^l--- (3.27) 

Thus, in this range of values of 5, we cannot associate with ^{t) a distribution of 
relaxation times (all the more since now Pjir) < 0). 

3.5. The noise spectral density 

Let us come back to the Mittag-Leffler memory kernel ^{t) as given by Eq. (3.1), 
and introduce its Fourier-Laplace transform defined as usual for real u by: 

/•oo 

7(a;)= / -i{t)e'''Ut. (3.28) 
Jo 

Inversely, we have: 

2 

7(t) = - / ^{uj) cos utdw. (3.29) 

71" Jo 

Eq. (3.29) displays the fact that it is enough to know 3?e7(a;) for a; > in order 
to get an integral representation for ^{t). Using the identity 7(0;) = 7(2; = —ioj)-, 
together with the expression (3.7) for ^{z), we get (assuming a; > 0): 

-/ ^ 5-1 sin(<^7r/2) + i [{utpf - cos(J7r/2)] 

7(0;) = 750; J a; > 0. (3.30) 

l + 2cos(57r/2)(a;to) +(wio) 

We have in particular: 

3?e 7(0;) = 75 sin(57r/2)a;'^~^ ^- ^, a; > 0. (3.31) 

l + 2cos(57r/2)(a;to)^ + (a;to) 

The noise spectral density (|F(a;)| ) defined for any real u as the Fourier 
transform of the correlation function {F{t)F{t')), reads, making use of the second 
fiuctuation-dissipation theorem (Eq. (1.2)): 

1 + 2cos(57r/2)(|a;|to) + (I^Ko) 
in accordance with the expression provided in [3]. 

If, instead of the Mittag-Leffler memory kernel, we use a memory kernel ^{t) 
about which we only know that it behaves at times t ^ to as an inverse power law 
of time (7(t » to) oc t~^), we usually write the noise spectral density as [1],[17]: 

{\F{u)f) = 2mkBT^s sin(57r/2)|a;|'^-' fc{\co\to), (3.33) 
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with /c a high-frequency cut-ofF function of typical width t^^. For 5 7^ 1, we 
have at hand a non-Ohmic dissipation model [17] (more precisely, sub-Ohmic for 
< 5 < 1, and super-Ohmic for 1 < 5 < 2). The spectral density of the coupling, 
defined for positive u by J{u>) = m^e^{u)), is given by: 

J{uj) = m^5sm{5n/2)uj^fc{ujto), uj > 0. (3.34) 

These results are in accordance with those obtained in [18], where it has been 
shown, in the framework of the Kac-Zwanzig model of a particle interacting 
through springs with a heat bath made of a collection of harmonic oscillators 
[19], how the distributions of frequencies, masses, and spring constants, determine 
the parameters of the generalized Langevin equation obeyed by the distinguished 
particle in the thermodynamic limit. The analysis carried out in [18] has been 
applied to the case of a memory kernel decreasing as an inverse power law of time. 

Let us come back to the full formula (3.32) for the noise spectral density. As 
put forward in [3], the choice of the Mittag-Leffler memory kernel allows us to 
make precise the cut-off function, which depends on the exponent S: 

fci\oj\to) = T 7^- (3.35) 

l + 2cos{S7r/2){\uj\tof + {\u;\tof' 

Note that, for 5=1, formula (3.35) for the cut-off function reduces to the Drude 
expression: 

fci\oj\to) = ^- (3.36) 

l + (|a;|to)' 

Interestingly, instead of being introduced by hand (usual choices being an exponen- 
tial or a Lorentzian cut-off function [17]) the 5-dependent cut-off function comes 
here into play in a natural way in the framework of the Mori-Zwanzig formalism 
(see Appendix B for more details). 

4. The particle's mean velocity 

The solution (f (t)) of the generalized Langevin equation (1.1) with 7(t) the Mittag- 
Leffler memory kernel as studied in Section 3 has first been obtained in [3] as the 
inverse Laplace integral: 

(„(,)) = „(« = 0)^/^^j^e"d.. (4.1) 
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with 7(2;) as given by Eq. (3.7). To get an explicit formula, we write: 



z + 'y{z) z 



(4.2) 



fc=0 



The expression in the r.h.s. of Eq. (4.2) can conveniently be rewritten as: 



E 

k=0 



(4.3) 



that is: 



E 

k=0 



k\ z^ V 



(-1) k\ /75 



5^ {z + t^'z^-^f^' to ^' {z + t-,'z^-^f^' 

(4.4) 

The result of the Laplace inversion takes the form of the sum of two infinite series 
of derivatives of generalized Mittag-Leffier functions (see for instance the formulas 
provided in [7]): 



./c=0 - "0 

doJ k 



.to) J 



-(;^)E 



(-1)^ /75\fc fcp(fc) 
k=o - ^0 



(4.5) 



The generalized Mittag-Leffier function Ea,j3{y) is defined through the series [6]- 
[7]: 



^-M!/) = Ef(:^ 



a > P>0, 



(4.6) 



and E^^p[y) represents its /cth-order derivative with respect to its argument y. 



(4.7) 



The main question we now want to address is that of the possible existence 
of a distribution of relaxation times associated with the particle's mean velocity. 
This discussion would be tricky to achieve on the basis of the full exact formula 
(4.5). This is why, from now on, we will use a simpler expression for {v(t)), valid 
for times Iq. 
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4.1. The average velocity at times to 

Let us thus take for 7(2;) its approximate expression (3.8) valid for \z\ <^ t^^, 
which amounts to use, in the time domain, the expression (3.5) for ^(t), vahd for 
times to. We then write accordingly (with obvious notations) [1]: 

(^(^))^ '"^^Zs^l N«^o\ (4.8) 

which gives: 

{v{t))'^v{t = ^)E2.s{-{y^6tf~\ t^to. (4.9) 

As shown in [3], this behaviour can also be retrieved from the general expression 
(4.5). Replace, in the two infinite series in the r.h.s. of Eq. (4.5), each general- 
ized Mittag-Leffler function by its asymptotic expression for t ^ to, using the 
asymptotic formula [7]: 

The leading order terms, which come from the second series in formula (4.5), can 
be resummed. The pertinent time scale is now Ug^ (the shorter time scale to being 
no more involved). We recover in this way formula (4.9). 

4.2. Connection with the fractional Langevin equation 

To shed some more light on the above results, let us remark that, using for ^{t) 
the approximate expression (3.5), amounts, in the case < 5 < 1, to writing the 
generalized Langevin equation (1.1) as a fractional Langevin equation, namely: 

mv(t) + m—^^[(t-t'y\{t')dt' = F{t), 0<S<1. (4.11) 

r(i - d) Jo 

In the range 1 < 5 < 2, we write instead, in order to avoid divergencies: 

mv{t)+m^s ''^l^^ +"^1X2^ ^ " fy'^'vit') dt' = F{t), 1< 5 < 2. 

(4.12) 

Accordingly, the average particle's velocity as given by Eq. (4.9) coincides as it 
should with the (averaged) solution of Eq. (4.11) (for < 5 < 1), or of Eq. (4.12) 
(for 1<5<2). 
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In addition, let us note that formulas (4.8) and (4.9) for the average particle's 
velocity have first been provided in [20] for < 5 < 1 as the force-free solution of a 
fractional Fokker-Planck equation, as well as in [21] for 1 < 5 < 2 as the force-free 
solution of a fractional Kramers equation, this latter equation describing both the 
velocity and coordinate of a particle exhibiting anomalous diffusion in an external 
force field (see also [9]). 

4.3. The noise spectral density 

Let us here quote for completeness the corresponding expressions for the noise 
spectral density and the spectral density of the coupling. In the angular frequency 
range |a;| <C t^^, formula (3.31) for 3?e7(a;) reduces to the expression used in [1], 
namely: 

SRe7(w) ~75sin(57r/2)|a;|^"\ (4.13) 
The corresponding noise spectral density (Eq. (3.33)) now reads: 

{\F{uj)f) ~ 2mkBT-fssm{67T/2)\ujf~'^. (4.14) 

As for the spectral density of the coupling (Eq. (3.34)), it reduces to: 

J{uj) c^m-fssm{57i/2)u^, u > 0. (4.15) 

Formulas (4.13), (4.14), and (4.15) are deduced from their counterparts valid for 
any u (respectively, formulas (3.31), (3.33), and (3.34)) by making fc{\uj\to) ~ 1 
(indeed the cut-off function does not play any role in the angular frequency range 

\oo\to < 1). 

4.4. Duality symmetry 

Let us compare formula (3.14) for 7(2;) (with ^{z) as given by Eq. (3.15)), with the 
approximate formula (4.8) for {v{z)). Both expressions have the same structure. 
They can even be identified, provided that one makes the change of indexes: 

5^2-5, (4.16) 

together with the change of time scales: 



(4.17) 
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This dynamical mapping between the Laplace transforms of the memory ker- 
nel, on the one hand, and of the particle's mean velocity, on the other hand, 
transforms the parameter region < 5 < 1 into the region 1 < 5 < 2 (and 
vice- versa), with important consequences about the existence of relaxation times 
distributions for either quantities. 

4.5. The relaxation times distribution associated with {v{t)) 

We will take advantage of the duality symmetry to discuss the existence of a 

relaxation times distribution associated with the particle's mean velocity as given 

by the approximate formula (4.9). In both cases < 5 < 1 and 1 < 5 < 2, we 

write: 

{v{t)) = Mt)) + {v2{t)), t>0, (4.18) 

where: 

{Mt))=v{t = 0):^ [ ^ e^' dz, (4.19) 

^TTZ iHa(£) Z + UJg ° Z^'^ 

and {v2{t)) is the contribution of the relevant poles of l/{z + z^~^), if any. 

In the case 1 < 5 < 2, no such poles exist, and the mean velocity relaxes 
towards zero without oscillating. We have: 

/•oo 

{v{t)) ^v{t^O) (r) e-*/^ dr, (4.20) 

^0 

with the relaxation times distribution (r) as given by: 

„ , , — sin^TT 1 1 / . X 

^ ^ (cosT) + 2 COS on + \u)sT) 

In the case < 5 < 1, the particle's mean velocity contains in addition an 
oscillating part: 



{v{t)) = v{t = 0) <j eU,,tcos[7r/(2-5)] cos[a;5tsin(^ 



2-6 



+ ^"p^,)(r)e-*/^dr|, 



(4.22) 

where P(i,)(t), although still given by Eq. (4.21), cannot obviously be interpreted 
as a relaxation times distribution. We easily check that, when S ^ 1, we have 
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— >■ 7i) and P(„)(t) ^ — (5(t — 7]7^)- Formula (4.22) thus consistently yields 
in this limit the standard behaviour corresponding to a non-retarded Langevin 
equation with friction coefficient 71: 

{v{t)) ^ v{t = 0) e-'^'\ (4.23) 

In the opposite limit S ^ 0, setting us — > uq, we have an undamped oscillation: 

{v{t)) v{t = 0) cosa;o^- (4.24) 

A physical explanation for the presence of an oscillating term in the r.h.s. of Eq. 
(4.22) can be provided in terms of the cage effect. As explained in [22], for small S 
the friction force induced by the environment is not just slowing down the particle, 
but also causes the particle a rattling motion. The generalized Langevin equation 
can be written, for < 5 < 1, as the fractional Langevin equation (4.11), that is: 

mv{t) + m-^^ [\uJs{t-t')]-'v{t')dt' = F{t). (4.25) 

r(i - 0) Jo 

In the limit 5 — > 0, Eq. (4.25) simply reads: 

mv{t) + mujl[x-x{t = 0)]=F{t). (4.26) 

This shows that a purely oscillating behaviour is expected in this limit. This os- 
cillation can be explained by the rattling motion of the particle of interest in the 
cage formed by the surrounding particles. 

4.6. Discussion 

Let us comment the above found results concerning the fact that it is not possible 
to define a relaxation times distribution for the particle's average velocity in the 
range < 5 < 1. 

Consider an assembly of independent particles obeying the usual simple Lange- 
vin equation: 

mv{t) + = F{t), (4.27) 

T 

where F{t) is the Langevin force and the relaxation time r is random. For any 
given T, the average velocity of an individual particle varies like v{t = 0) exp(— t/r). 
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An obvious question is whether it exists a relaxation times distribution function 
P{v) ("T") such that the average velocity as defined by: 

/"OO 

{v{t)) = v{t = 0) / P(„) (r) e-*/- dr, (4.28) 
effectively obeys the (averaged) fractional Langevin equation, namely: 

rn{v{t)) +m—^^ / {t - t'y\v{t')) dt' ^ 0, < 5 < 1. (4.29) 

r(i - d) Jo 

It ensues from the results of Subsection 4.5 that such a distribution function cannot 
be defined. In other words, it does not exist a relaxation times distribution which 
would map the simple Langevin equation onto the fractional one. 

5. The scattering function 

We continue in this section, and in the following of this paper as well, to study the 
various functions of time which enter into play at times t ^ to only. This latter 
time scale will thus no more appear in what follows. 

5.1. The mean- squared displacement 

Let us introduce the particle displacement x{t) ~ jQv{t')dt'. The two-time dis- 
placement correlation function can be obtained from the velocity correlation func- 
tion as the double integral: 

{x{t)x{t'))^ [ [ {v{h)v{t2)) dhdt2. (5.1) 

JQ Jo 

On the basis of the first fluctuation-dissipation theorem (see formula (B.6)) to- 
gether with formula (4.9) for the particle's mean velocity, the velocity correlation 
function is given by [1] , [9] : 

{v{t)v{t')) = ^E2-s [-Mt - t'\f-'] . (5.2) 
m 

To obtain {x{t)x{t')), we have to use the integration formulas: 

/ E2-s[-{^5ti) ]dh^tE2-5,2[-ioostf~^], (5.3) 
Jo 

and: ^ ^ 

/ dh [ ' E2-s[-{^5t2y~']dt2 ^t^ E2-s,3[-{ujstf~']. (5.4) 
Jo Jo 
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Making t' = t in Eq. (5.1), gives the expression for the mean-squared displacement 
(Ax^t)) = {x\t)) [1],[9],[21]: 

{Ax\t)) ^2^e E2-s,3 . (5.5) 

In the case < 5 < 1, the expression for {Ax'^{t)) contains an oscillating 
part. The general formula being a bit lengthy, we shall only quote its two limiting 
forms for 5 — > 1 and 5 — > 0. First, we can check that, when 5 — > 1, we recover the 
standard Brownian motion result: 



In the opposite limit 5 — )■ 0, we get the undamped oscillator formula: 



(Ax^it)) = 2^ (1 - cosujot). (5.7) 



In the case 1 < 5 < 2 in which the motion is superdiffusive, we can write the 
mean-squared displacement in terms of the relaxation times distribution function 
P{v) (t) as: 

7 rji /»00 

^AxHt))=2^ / P^,^{T)[t-T{l-e-'/n]TdT. (5.8) 
™ Jo 

5.2. Large time behaviour of the mean-squared displacement: the non- 
inert ial regime 

We mean here by "large times" the times for which we have uigt ^ 1. The general- 
ized Mittag-Leffler function in Eq. (5.5) can then be replaced by its leading order 
term as given by the asymptotic formula (4.10). We get [1]: 

Formula (5.9) is valid for any 5 {0 < d < 2). Note that, when u)st S> 1, the 
oscillating contribution which exists at shorter times in the range < 5 < 1 does 
not appear any more. 

Interestingly, the result (5.9) is fully consistent with the expression for the 
velocity correlation function obtained by considering, instead of the full generalized 
Langevin equation (1.1), the corresponding noninertial equation of motion, i.e.: 



m 



[ -f{t-t')v{t')dt' ^ F{t). (5.10) 
Jo 
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Applying Laplace transformation to Eq. (5.10), we get: 

m^{z)v{z) = F{z), (5.11) 

hence the expression for the double Laplace transform of the velocity correlation 
function: 

irm) = (5.12) 
From the second fluctuation-dissipation theorem, we have [1] : 

{F{z)F{z')) = mfegT ^^^^^y^ (5.13) 

z + z 

and thus: 

{v{z)v{^')) = ^^(^ + T7^)- (5.14) 

^ ^ ^ ^ ^' m z + z' \'y{z) 'y{z')J ^ ' 

Taking for 7(2;) its expression (3.9) valid for \z\ <^ ^, we get: 

= ^||^(^'-^ + ^"-*). (5.15) 

Neglecting the inertial term thus yields for the velocity correlation function the 
expression: 

^'^'^'^ = W^)' ^ ^ ^ 

This latter formula is identical to the one which can be deduced by derivating 
twice with respect to time the large time expression (5.9) for (Ax^(t)). This fact is 
in accordance with the general view according to which considering the noninertial 
equation of motion is tantamount to studying the motion of the particle over a 
coarse-grained time scale (here, u^^). 

5.3. The scattering function 

Let us introduce the scattering function Fs{q.,t), as defined by [23]: 

Fs{q,t) = J Gs{x,t)e'^^dx, (5.17) 

with Gs{x,t) the diffusion profile (if the present one-particle model is considered 
as picturing the diffusion of an ensemble of non-interacting particles, the function 
Fs{q, t) describes the relaxation of a density fiuctuation of wavevector q). Since the 
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noise has been assumed to be Gaussian, the diffusion profile is a Gaussian function 
of width {Ax^{t)): 

„2 



Gs{x,t) = [27r{Ax^{t))] ^/^exp 



X 



2{Ax'^(t)). 

Performing the Fourier integration in Eq. (5.17), we get: 

,{Ax\t))- 



(5.18) 



Fs{q,t) = exp 



-Q 



(5.19) 



Using Eq. (5.5), which gives the expression for the mean-squared displacement 
in terms of a generalized Mittag-Leffler function, would lead to a fairly involved 
formula for Fs{q,t). Let us instead simply use for {Ax^{t)) its expression (5.9) 
valid for ust ^ 1 (non-inertial regime). The scattering function is then simply a 
modified exponential function [24]: 



Fs{q,t) ~ exp 



2 fcgT {ust) 
^ ma;|r(5+l) 



u^t » 1. 



(5.20) 



5.4. The relaxation times distribution associated with Fs{q,t) 

In the range < 5 < 1 which corresponds to subdiffusion, the scattering 
function as given by formula (5.20) is a stretched exponential, whereas, in the range 
1 < 5 < 2 in which the motion is superdiffusive, it is a compressed exponential. 
Accordingly, a relaxation times distribution can be associated with the scattering 
function in the range < 5 < 1, but not in the range 1 < S < 2. 
Introducing the time scale tq as defined by: 



1 



rmoi 



(5.21) 



we can rewrite Eq. (5.20) as: 

Fs{q,t) ~ exp 



Lost > 1. 



(5.22) 



r(5 + l)J 

Formula (5.22) for Fs{q,t) has the same form as the approximate formula (3.3) 
for 7(t) valid for t <^to- Both formulas can be identified (up to a proportionality 
factor), provided that one makes the change of time scales: 



to i — > tq, 



(5.23) 
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the parameter 5 remaining unchanged. 

In the range < S < 1, introducing the rates distribution fiF^{k), we write: 



Fs{q,t) 



lJiF,{k)e dk 



that is, in terms of the relaxation times distribution Pp^{r): 

Fs{q,t)^ / Pf^{t) e-'/^ dr. 
Jo 

The relaxation rates distribution jJiPsi^) in Eq. (5.24), as given by: 



2771 



exp 



Br 



ns+i) 



dz, 



identifies with the one-sided stable Levy law ^{k) as defined by [25]: 



(5.24) 



(5.25) 



(5.26) 



(5.27) 



provided that one sets C = /V[5 + 1). This property deserves to be mentioned 
since it relates the fractional Langevin equation to the generalized central limit 
theorem [25]. 

Carrying out a discussion similar to the one in Appendix A, we write the 
relaxation rates distribution /iF^ik) as: 



I^F, (k) 



1 r 



exp 



sm 



r(5 + i)J Lr((5 + i)J 



e ''^"''dx. 



As for the relaxation times distribution Pf^{t), it reads: 



It 

PFAr) = -:^,j^ exp 
For T <^ tg, we simply have: 



cos Stt 



''x^ sin 5tt 



sm 



r(5 + i) 



e-^'-^/^dx. 



(5.28) 



(5.29) 



(5.30) 



In the range 1 < 5 < 2, Fg{q^t) is a compressed exponential with which no 
relaxation times distribution can be associated. Finally, let us note for complete- 
ness that, in the limiting case 5 = 1, we just recover an exponentially decaying 
scattering function: 



Fs{q, t) ~ exp(-t/t<,), 71^ > 1, 



with: 



1 2kBT 



t 



m7i 



(5.31) 
(5.32) 



6. Conclusion 
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In this paper, we studied the dynamics of an anomalously diffusing particle whose 
motion is described by a generalized Langevin equation and whose mean-squared 
displacement grows at large time like {Ax'^{t)) ~ t^, with < 5 < 1 (subdiffu- 
sion) or 1 < 5 < 2 (superdiffusion). We successively focused our interest on three 
time decaying quantities, namely, the memory kernel of the generalized Langevin 
equation, the particle's mean velocity, and the scattering function. 

Using for the memory kernel ^{t) a Mittag-Leffler function of index 6 and 
characteristic time scale to, we showed that its decay can be associated with a 
relaxation times distribution in the range < 5 < 1. We were able to make 
precise this distribution at any time t, as well as to check its consistency with the 
relaxation times distributions associated with the limiting forms of the Mittag- 
Leffier function for t to {a stretched exponential) and for t ^ to (an inverse 
power law of time). Interestingly, the distribution of small relaxation times r ^ to 
was shown to exhibit an integrable divergence oc t^~^. As stated in [24], such a 
divergence implies a profusion of short times scales in the Mittag-Leffler decay of 
index 5 in the range < 5 < 1. As for the distribution of large relaxation times, it 
was shown to decay like an inverse power law oc , which displays the fact that 
there exist also abundant long time scales in the Mittag-Leffler decay. 

We then turned to the particle's mean velocity {v{t)). Restricting the study to 
times t^ to, we showed that {v{t)) is given by a Mittag-Leffler function of index 
2 — 5 and characteristic time scale » to, and thus that it can be associated 
with a relaxation times distribution in the range 1 < 5 < 2. More precisely, the 
dynamical mapping which allows us to deduce the behaviour of {v{t)) from that of 
7(t) was shown to correspond to the change of indexes 5 i — > 2 — S, together with 
the change of time scales to < — > ^s^^- The divergence of the distribution of small 
relaxation times of 7(t), all the more pronounced as 6 approaches zero, induces 
for < t < to a,n efficient decrease of the dissipative kernel, resulting in turn in 
a weak damping of the particle's mean velocity, which is actually described by an 
oscillating function in the range < d < 1. Conversely, for 1 < 6 < 2, the decrease 
of 7(t) for < t < to is relatively weak, which in turn implies an efficient decrease 
of the particle's mean velocity. 

Finally, we turned our interest on the scattering function Fs{q,t), whose be- 
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haviour is directly linked to that of the mean-squared displacement. It results from 
the very definition of the scattering function that this quantity can be associated 
with a relaxation time distribution in the subdiffusive range < S < 1, but not in 
the superdiffusive range 1 < 5 < 2. 
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Appendix A 

Let us check the consistency, in the case < d < 1, between the relaxation times 
distributions obtained for the Mittag-Leffler function Es[—{t/to)^], on the one 
hand, and for the stretched exponential function and the algebraically decreasing 
function respectively describing its short time and long time behaviours, on the 
other hand. 



A.l. The relaxation times distribution associated with the stretched 
exponential exp[-|^^^] (0 < 5 < 1) 

Let us here discuss the existence of a distribution of relaxation times for the mod- 
ified exponential function exp[— which pictures the short time behaviour 
of the Mittag-Leffler function Es[-{t/toY]. 

The stretched exponential exp[— r(g+i) ] with < (5 < 1 is completely mono- 
tonic (it is often referred to as the Kohlrausch- Williams- Watt (KWW) relaxation 

(t/t 

law [26]), but this is not the case for the compressed exponential exp[— r/j+i) ] with 
1 < 5 < 2. It is therefore solely in the range < 5 < 1 that a relaxation times dis- 
tribution, denoted by -Pkww(^): can be associated with the modified exponential 



exp 







I r{s + i)\ 





A*KWw(^)e dk 



(A.l) 



In Eq. (A.l), the relaxation rates distribution //kww(^) ^ given by: 
/.Kww(/c) = ^/ exp[-— ^Je dz, 



(A.2) 



identifies with the one-sided stable Levy law ^{k) with C = /T{S + 1). We 
have [4], [26]: 



1 

//Kww(/c) ^ - to exp 

71" Jo 



COS 5n 



[x^ sin Stt 



sm 



T{5+1) 



mi 



e ^^""^ dx. 



(A.3) 



hence the associated relaxation times distribution: 



KWW 



(r) 



1 to 



exp 



x^ cos (Jtt 



sm 



r(5 + i)J Lr(5+i)J 



x^ sin ^TT 
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Using a Tauberian theorem [4], it is easy to deduce the behaviour of Pkww{'t) 
for T <^to from the smaU-a; behaviour of the function exp[-^^^p^] sin[^^|||y^]. 
As a result, we get: 

PKww{r) ^ r(rr) ' T < to- (A.5) 

TT to \to/ 

This result is fully consistent with the form taken by P-yir) (i.e. the distribution of 
relaxation times associated with the Mittag-LefHer function Es[—{t/to)^]) in the 
same range of values of r (see Eq. (3.25)). Formulas (3.25) and (A.5) display an 
integrable divergence of the distribution of small relaxation times, which implies a 
profusion of short time scales in both the Mittag-Leffler function and the associated 
stretched exponential [24]. 

A. 2. The relaxation times distribution associated with the algebraically 
decaying function ^^^^^ (0 < 6 < 1) 

In the same range of values of 6, consider the algebraic function ^i-^) which 
pictures the long time behaviour of the Mittag-Leffler function Eg[—{t/to) ]. In- 
terestingly, this algebraic function can be written as: 



Eq. (A. 6) displays the fact that the associated relaxation times distribution is 

sin^TT 1 / T \ 



TT T Vt 



(A.7) 



a result consistent with the form taken by the function P^(t) as given by Eq. (3.23) 
for large relaxation times (i.e., for r ^ to). 



Appendix B 
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B.l. The relation between ^{z) and (l){z) 

Let us here come back in more detail to the relation (3.14) between 7(2;) and (l>{z), 
rewritten below for clarity: 

= (B.1) 

z + (p[z) 

together with the generalized Langevin equation (1.1), rewritten in Laplace trans- 
forms notations for convenience: 

zv{z) - v{t = 0) + j{z)v{z) = —Fi{z), (B.2) 

m 

and the Langevin- like equation (3.10) for Fi{t), also rewritten in Laplace trans- 
forms notations: 

zFi{z) - Fi(t = 0) + ^{z)Fi{z) = F2{z). (B.3) 

Following Mori's procedure [11], we write the particle's velocity as a sum of 
two orthogonal components, namely, its projective and vertical components with 
respect to the v axis (the symbol v denotes the value of v{t) at time t = 0): 

.(t) = MM.+ /*MiA(t-t')*'. (B.4) 

[v ) Jo ) 

or, in Laplace transforms notations: 

= ^^s^ ^ + -hi-)- (B-5) 

Identifying Eqs. (B.2) and (B.5), gives: 

11/"°° . , 

{v{t)v)e''Ut, h{z)^Fi{z), (B.6) 



z + 7(2) {v^) Jo 



i.e. the first fluctuation-dissipation theorem. In the same way, still following Mori's 
procedure, we write the Langevin force as a sum of its projective and vertical 
components with respect to the Fi axis: 

F,{t)-^^F, + -^^-f,{t-t)dt, (B.7) 
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or, in Laplace transforms notations: 



F,{z) = ^» ^ F, + ^° -f2{z). (B.8) 



(B.9) 



Identifying Eqs. (B.3) and (B.8), gives: 

11/"°° . . 

ITT = 7P2Y / {Fi{t)Fi)e-^' dt, h{z) = F2{z) 

z + (j){z) {J^i) Jo 

Now, the second fluctuation-dissipation theorem reads: 

{Fr{t)Fr) = mkBT-i{t), (B.IO) 

or, in Laplace transforms notations: 

/ {Fi{t)Fi)e-'Ut^mkBT^{z). (B.ll) 
Jo 

Comparing Eqs. (B.9) and (B.ll), gives: 

^{z)= \ ■ (B.12) 

Then, making use of formula (B.IO) at time t = 0, we recover formula (B.l). 
B.2. Derivation of the noise cut-oflF function 

Interestingly, in the framework of the Mori-Zwanzig formalism, the noise cut-off 
function is not introduced a priori for convenience purposes, but can be consistently 

derived once the second-order memory kernel (f){t) is known. 

Applying standard harmonic analysis to Eq. (3.10), we get, introducing the 
Fourier-Laplace transform (l){u)) of (pit): 



{\F,{uf) = ^'-^^^7' o - (B.13) 
\—i(jO -\- (j>{(jo)\ 

On the other hand, the second- fluctuation theorem (1.2) reads, using Fourier trans- 
forms: 

{\Fi{uj)\^) = 2mkBTne%u). (B.14) 



1-^2 (^0 



Since we have: 
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Eq. (B.14) reads: 



^0 + 



From formulas (B.13) and (B.16), we deduce the spectral density of the second- 
order random force: 

{\F2{oo)f) = 2mkBT^^e4>{oo), (B.17) 

*o 

that is, making use of the identity = (f){z = —iio), together with the expression 
(3.15) for 4){z): 

(|F2(a;)|^) = 2mA;BT75sin(57r/2)to^Vl^"'^- (B.18) 

Let us now come back to the spectral density (|Fi(a;)|^) as given by Eq. (B.16). 
We have: 

{\F,iu)\') = 2mksT:^^^M^ , ^, (B.19) 

an expression, which, once the expression for ^e(f){u>) is made explicit, identifies 
with formula (3.33) of the main text, that is: 

(|Fi(a;)|') = 2mkBT^s sm{57r /2)\uj\^-^ fc{\uj\to). (B.20) 

The shorter time scale to intervenes only in the cut-off function. This latter verifies: 

/.(|.|*„) = "''1""°' (B.21) 
\—iu) + 0(a;)| 
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